In search of quark gluon plasma in nuclear collisions 
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Abstract 

At high temperatures and densities the nu- 
clear matter undergoes a phase transition to a 
new state of matter called quark gluon plasma 
(QGP). This new state of matter which ex- 
isted in the universe after a few microsecond of 
the big bang can be created in the laboratory 
by colliding two nuclei at relativistic energies. 
In this presentation we will discuss how the 
the properties of QGP can be extracted by an- 
alyzing the spectra of photons, dileptons and 
heavy flavours produced in nuclear collisions 
at Relativistic Heavy Ion Collider (RHIC) and 
Large Hadron Collider (LHC) energies. 



I. INTRODUCTION 

The theory of strong interaction - Quantum 
Chromodynamics (QCD) has a unique feature 
- it possess the property of asymptotic freedom 
which implies that at very high temperatures 
and/or densities nuclear matter will convert 
to a deconfined state of quarks and gluons [1]. 
Recent lattice QCD based calculations 2] in- 
dicate that the value of the temperature for 
the nuclear matter to QGP transition ^ 175 
MeV. It is expected that such high tempera- 
ture can be achieved in the laboratory by col- 
liding nuclei at RHIC and LHC energies. 

A high multiplicity system of deconfined 
quarks and gluons with power law type of mo- 
mentum distributions can created just after 
the nuclear collisions at high energies. Interac- 
tions among these constituents may alter the 
momentum distribution of quarks and gluons 
from a power law to an exponential one - re- 
sulting in a thermalized state of quarks and 
gluons with initial temperature, T^. This ther- 
malized system with high internal pressure ex- 
pands very fast as a consequence it cools and 



reverts to hadronic matter at a temperature, 
Tc ^ 175 MeV. The hadrons formed after the 
hadronization of quarks may maintain thermal 
equilibrium among themselves until the ex- 
panding system becomes too dilute to support 
collectivity at a temperature, Tp(~ 120 MeV) 
called freeze-out temperature from where the 
hadrons fly freely from the interaction zone to 
the detector. 

The electromagnetic (EM) probes [3J (see 
"41-16!] for review) i.e. real photons and dilep- 
tons can be used to follow the evolution of the 
system from the pristine partonic stage to the 
final hadronic stage through an intermediary 
phase transition or cross over. In the state 
of QGP some of the symmetries of the phys- 
ical vacuum may either be restored or bro- 
ken - albeit transiently. The electromagnetic 
probes, especially the lepton pairs can be used 
very effectively to investigate whether these 
symmetries in the system are restored/broken 
at any stage of the evolving matter. Results 
from theoretical calculations will be shown in 
the presentation to demonstrate this aspect 
of the electromagnetically interacting probes. 
We will demonstrate that lepton pairs can be 
used very effectively to probe the collective 
motion (radial and elliptic) of the system. 

The other promising probe of the QGP that 
will be discussed here - is the depletion of 
the transverse momentum spectra of energetic 
quarks (and gluons) in QGP. The magnitude 
of the depletion can be used to estimate the 
transport coefficients of QGP which is turn 
can be used to understand the fluidity of the 
matter. 

The transport coefficients of QGP and 
hot hadrons calculated by using perturbative 
QCD and effective field theory respectively 
have been applied to evaluate the nuclear sup- 
pression {Raa) of heavy fiavours. Theoreti- 
cal results on Raa will be compared with the 



experimental data available from RHIC and 
LHC energies. The azimuthal asymmetry of 
the system estimated through the single lep- 
tons originating from the decays of open heavy 
flavours produced from the fragmentation of 
heavy quarks will also be discussed. 



II. THE ELECTROMAGNETIC 
PROBES 

The dilepton production per unit four- 
volume from a thermal medium produced in 
heavy ion collisions is well known to be given 
by: 
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where the factor L{AP) = {l + 2mf/A'P) (1- 
'imf/M'^y/^ is of the order of unity for elec- 
trons, M(= \/p^) being the invariant mass of 
the pair and the hadronic tensor PV^y is de- 
fined by 
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where Jf^{x) is the electromagnetic current 
and (.) indicates ensemble average. For a de- 
confined thermal medium such as the QGP, 
Eq. ([T]) leads to the standard rate for lepton 
pair productions from qq annihilation at low- 
est order. 

The production of low mass dileptons from 
the decays of light vector mesons in the 
hadronic matter can be obtained as (see 7] 
for details): 
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where /be is the thermal distributions for 
bosons, p^y{qo,q) is the spectral function of 
the vector meson R{= PtUJ, </>) in the medium, 

^R = ^r'^'r^ ^^R is the mass of R and Fr is 
related to the decay of R to lepton pairs. 



The interaction of the vector mesons with 
the hadrons in the thermal bath will shift the 
location of both the pole and the branch cuts 
of the spectral function - resulting in mass 
modification or broadening - which can be 
detected through the dilepton measurements 
and may be connected with the restoration of 
chiral symmetry in the thermal bath. In the 
present work the interaction of p with thermal 
7r,aj, ai, hi [7|, |8| and nucleons [2] have been 
considered to evaluate the in-medium spectral 
function of p. The finite temperature width of 
the Lu spectral function has been taken from 

To evaluate the dilepton yield from a dy- 
namically evolving system produced in heavy 
ion collisions (HIC) one needs to integrate the 
fixed temperature production rate given by 
Eq.[3]over the space time evolution of the sys- 
tem - from the initial QGP phase to the fi- 
nal hadronic freeze-out state through a phase 
transition in the intermediate stage. We as- 
sume that the matter is formed in QGP phase 
with zero net baryon density at temperature 
Ti in HIC. Ideal relativistic hydrodynamics 
with boost invariance [ll| has been applied to 
study the evolution of the system. 

The EoS required to close the hydrody- 
namic equations is constructed by taking re- 
sults from lattice QCD for high T [2| and 
hadron resonance gas comprising of all the 
hadronic resonances up to mass of 2.5 GeV 
[i3,[l3] for lower T. The system is assumed to 
get out of chemical equilibrium at T = Teh — 
170 MeV [14]. The kinetic freeze-out tempera- 
ture Tp = 120 MeV fixed from the pT spectra 
of the produced hadrons. 



A. Invariant mass spectra of lepton pairs 

The M distribution of lepton pairs origi- 
nating from quark matter (QM) and hadronic 
matter (HM) with and without medium ef- 
fects on the spectral functions of p and uj 
are displayed in Fig. [1] We observe that for 
M > Mcj, the QM contributions dominate. 
For Mp < M < M^ the HM shines brighter 
than QM. For M < Mp, the HM (solid fine) 



over shines the QM due to the enhanced con- 
tributions primarily from the medium induced 
broadening of p spectral function. However, 
the contributions from QM and HM become 
comparable in this M region if the medium 
effects on p spectral function is ignored (dot- 
ted line). Therefore, the results depicted in 
Fig.[T]indicate that a suitable choice of M win- 
dow will enable us to unravel the contributions 
from a particular phase (QM or HM). An ap- 
propriate choice of M window will also allow 
us to extract the medium induced effects. 

To further quantify these points we evaluate 
the following |15J |: 
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where the M integration in both the numera- 
tor and denominator are performed for selec- 
tive windows from Mi to M2 with mean M de- 
fined as (M) = {Mi + M2)/2. While in the de- 
nominator the integration is done over the en- 
tire lifetime, the prime in J in the numerator 
indicates that the r integration in the numer- 
ator is done from ri = r^ to T2 = Ti + At with 
incremental At until T2 attains the life time of 
the system. In the inset of Fig. [T]i^ is plotted 
against Tav{= (n -I- T2)/2). The results sub- 
stantiate that pairs with high (M) ^ 2.5 GeV 
originate from early time {r^^y < 3 fm/c, QGP 
phase) and pairs with (M) ~ 0.77 GeV mostly 
emanate from late hadronic phase (ra^ > 4 
fm/c) . The change in the properties of p due 
to its interaction with thermal hadrons in the 
bath is also visible through F evaluated for 
(M) ~ 0.3 GeV with and without medium ef- 
fects. 



B. Elliptic flow of lepton pairs 

The elliptic flow of dilepton, V2{pt, M) can 
be defined as: 



V2 = 
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FIG. 1: Invariant mass distribution of dileptons 
from hadronic matter (HM) for modified and un- 
modified p meson. 
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FIG. 2: [a] and [b] indicate elliptic flow of lepton 
pairs as a function of pr for various M windows. 
[c] displays the effect of the broadening of p spec- 
tral function on the elliptic flow for (M) = 300 
MeV. [d] shows the variation of 7? (see text) with 
PT for (M) = 0.3 GeV, 0.77 GeV and 2.5 GeV. 



where the ^ stands for summation over QM 
and HM phases. 

Fig-[2l(N and [b]) show the differential el- 
liptic flow, V2{jpt) of dileptons arising from 
various (M) domains. We observe that for 
(Af) =2.5 GeV V2 is small for the entire px 
range because these pairs arise from the early 
epoch (see inset of Fig. [T|) when the flow is 
not developed entirely. However, the V2 is 
large for {M) = 0.77 GeV as these pairs orig- 
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FIG. 3: Variation of dilepton elliptic flow as 
function of (A/) for QM, HM (with and with- 
out medium effects) and for the entire evolution. 
The inset shows the variation of momentum space 
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pict the variation of R with px for (M) = 0.3 
GeV (solid circle) 0.77 GeV (solid line) and 
2.5 GeV (open circle), the quantity R is de- 
fined as _R = ^2 /{v2 + w^^^) where Wj is 
the elliptic flow of the phase i{= QM + HM. 
The results clearly illustrate that V2 of lepton 
pairs in the large (M) domain originate from 
QM. 

Fig. [3] shows pt integrated elliptic flow, 
V2{{M)) evaluated for different (M) windows 
defined above. The V2 (which is proportional 
to momentum space anisotropy, Cp) of QM 
is small because the pressure gradient is not 
fully developed in the QGP phase as evident 
from the inset plot of Cp with r. The hadronic 
phase V2 has a peak around p pole indicat- 
ing large flow at late times. For (M) > m^ 
the V2 obtained from the combined phases ap- 
proach the value corresponding to the V2 for 
QGP. Therefore, measurement of V2 for large 
(M) will bring information of the properties 
of the QGP. It is important to note that the 
Pt integrated V2{{M)) of lepton pairs with 
(M) ~ mTr,mK is close to the hadronic wj 
and V2 if the thermal effects on p properties 
are included. Exclusion of medium effects give 
lower V2 for lepton pairs compared to hadrons. 
We also observe that the variation of V2{{M)) 
with (M) has a structure similar to dN/dM 
vs M. As indicated by Eq. [1] we can write 
V2{{M)) - E4x/i- The structure of dA^/dAf 
is reflected in V2{{M)) through fi. 



FIG. 4: The dilepton yield plotted against Mt — 
Mav for different M windows for LHC initial con- 
dition. 



C. Radial flow of dileptons 



inate predominantly from the late hadronic 
phase when the flow is fully developed. It 
is also interesting to note that the medium 
induced enhancement of p spectral function 
provides a visible modification in V2 for dilep- 
tons below p peak (Fig. [2] [c]). The medium- 
induced effects lead to an enhancement of V2 
of lepton pairs which is culminating from the 
'extra' interaction (absent when a vacuum p 
is considered) of the p with other thermal 
hadrons in the bath. In Fig. [5] [d] we de- 



The transverse mass distributions of the lep- 
ton pairs at LHC is displayed in Fig. |31 The 
variation of inverse slope (deduced from the 
from the transverse mass distributions of lep- 
ton pairs, Fig. S]) with {M) for LHC is de- 
picted in Fig. [S] The radial flow in the sys- 
tem is responsible for the rise and fall of T^fj 
with {M) (solid line) in the mass region (0.5 < 
M(GeV)< 1.3), for wt = (dashed line) a 
completely different behaviour is obtained. 
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FIG. 5: Teff for difTerent values of the M-bins for 
LHC conditions. The dashed line is obtained by 
setting vt = 0- 



D. Radial flow from HBT interferometry 
of lepton pairs 

It was shown in Ref. [Ig that the varia- 
tion of HBT radii {Rside and Rout) extracted 
from the correlation of dilepton pairs with 
(M) can used to extract collective proper- 
ties of the evolving QGP. While the radius 
{Rside) corresponding to qside is closely re- 
lated to the transverse size of the system and 
considerably affected by the collectivity, the 
radius {Rout) corresponding to qout measures 
both the transverse size and duration of par- 
ticle emission. The extracted Rside and Rout 
for different {M) are shown in Fig. [6] The 
Rside shows non-monotonic dependence on M, 
starting from a value close to QGP value 
(indicated by the dashed line) it drops with 
increase in M finally again approaching the 
QGP value for (M) > m^. It can be shown 

that Rside ^ 1/(1 + EcoUective / Ethermal) ■ In 

the absence of radial flow, -Rside is indepen- 
dent of Qside- With the radial expansion of 
the system a rarefaction wave moves toward 
the center of the cylindrical geometry as a con- 
sequence the radial size of the emission zone 
decreases with time. Therefore, the size of 
the emission zone is larger at early times and 
smaller at late time. The high (M) regions are 
dominated by the early partonic phase where 
the collective flow has not been developed fully 
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FIG. 6: -Rside and Rout as a function of {M). The 
dashed, dotted and the solid line (with asterisk) 
indicate the HBT radii for the QGP, hadronic and 
total dilepton contributions from all the phases 
respectively. The solid circles are obtained by 
switching off the contributions from p and uj. 



i. e. the ratio of collective to thermal energy is 
small hence show larger Rside for the source. 
In contrast, the lepton pairs with M ~ rup are 
emitted from the late hadronic phase where 
the size of the emission zone is smaller due to 
larger collective flow giving rise to a smaller 
Rside- The ratio of collective to thermal en- 
ergy for such cases is quite large, which is re- 
flected as a dip in the variation of Rside with 
{M) around the p-mass region (Fig. [6] upper 
panel). Thus the variation of -Rside with M 
can be used as an efficient tool to measure the 
collectivity in various phases of matter. The 
dip in Rside at (Af) ^ rup is due to the contri- 
bution dominantly from the hadronic phase. 
The dip, in fact vanishes if the contributions 
from p and uj is switched off (circle in Fig. ^ . 
We observe that by keeping the p and to con- 
tributions and setting radial velocity, Vr — 0, 
the dip in -Rside vanishes, confirming the fact 
that the dip is caused by the radial flow of 
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FIG. 7: Raa as a function of pT for _D and B 
mesons at LHC. Experimental data taken from 



the hadronic matter. Therefore, the value of 
Rside at (M ) ^ TTip may be used to estimate 
the average Vr in the hadronic phase. 

The Rout probes both the transverse dimen- 
sion and the duration of emission as a conse- 
quence unhke Rside it does not remain con- 
stant even in the absence of radial flow and 
its variation with M is complicated. The large 
M regions are populated by lepton pairs from 
early partonic phase where the effect of flow is 
small and the duration of emission is also small 
- resulting in smaller values of Rout- For lep- 
ton pair from AI '^ nip the flow is large which 
could have resulted in a dip as in Rside in this 
M region. However, Rout probes the dura- 
tion of emission too which is large for hadronic 
phase. The larger duration compensates the 
reduction of Rout due to flow in the hadronic 
phase resulting in a bump in Rout for M ~ nip 
(Fig. |6] lower panel). Both Rside and Rout ap- 
proach QGP values for {M) - 2.5 GeV im- 
plying dominant contributions from partonic 
phase. 



duced in proton -I- proton (pp) collisions has 
been considered as a signature of QGP for- 
mation. The two main processes which cause 
the depletion are (i) the elastic collisions and 
(ii) the radiative loss or the inelastic collisions 
of the high energy partons with the quarks, 
anti-quarks and gluons in the thermal bath. 



In the present work we focus on the en- 
ergy loss of heavy quarks in QGP in deduc- 
ing the properties of the medium. Because (i) 
the abundance of charm and bottom quarks 
in the partonic plasma for the expected range 
of temperature to be attained in the exper- 
iments is small, consequently the bulk prop- 
erties of the plasma is not decided by them 
and (ii) they produce early and therefore, can 
witness the entire evolution history. Hence 
heavy quarks may act as an efficient probe 
for the diagnosis of QGP. The depletion of 
heavy quarks in QGP has gained importance 
recently in view of the measured nuclear sup- 
pression in the pt spectra of non-photonic sin- 
gle electrons [iTlflst. 



We assume here that the light quarks and 
gluons thermalize before heavy quarks. The 
charm and bottom quarks execute Brownian 
motion [19] (see references therein) in the heat 
bath of QGP. Therefore, the interaction of the 
heavy quarks with QGP may be treated as 
the interactions between equilibrium and non- 
equilibrium degrees of freedom. The Fokker- 
Planck (FP) equation provide an appropri- 
ate framework for the evolution of the heavy 
quark in the expanding QGP heat bath which 
can be written as [19| : 



III. SUPPRESSION OF HEAVY 
FLAVOURS IN QGP 

The depletion of hadrons with high trans- 
verse momentum {pt) produced in Nucleus + 
Nucleus collisions with respect to those pro- 
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where the momentum, p = (Pt^Pz) — 
{px,Py,Pz), 7 is the drag coefficient and D 
is the diffusion coefficient. We numerically 
solve Eq.[S] [201] with the boundary conditions: 
f{px,Py,t) — >■ for Px,Py —5- 00 and the ini- 
tial (at time t = Ti) momentum distribution 
of charm and bottom quarks are taken MNR 
code [m. 

The system under study has two compo- 
nents. The equilibrium component, the QGP 
comprising of the light quarks and the glu- 
ons. The non-equilibrium component, the 
heavy quarks produced due to the collision 
of partons of the colliding nuclei has momen- 
tum distribution determined by the perturba- 
tive QCD (pQCD), which evolves due to their 
interaction with the expanding QGP back- 
ground. The evolution of the heavy quark 
momentum distribution is governed by the FP 
equation. The interaction of the heavy quarks 
with the QGP is contained in the drag and dif- 
fusion coefficients. The drag and diffusion co- 
efficients are provided as inputs, which are, in 
general, dependent on both temperature and 
momentum. The evolution of the tempera- 
ture of the background QGP with time is gov- 
erned by relativistic hydrodynamics. The so- 
lution of the FP equation at the (phase) tran- 
sition point for the charm and bottom quarks 
gives the (quenched) momentum distribution 
of hadrons {B and D mesons) through frag- 
mentation. The fragmentation of the initial 
momentum distribution of the heavy quarks 
results in the unquenched momentum distri- 
bution of the B and D mesons. The ratio of 
the quenched to the unquenched pt distribu- 
tion is the nuclear suppression factor which 



is experimentally measured. The quenching 
is due to the dragging of the heavy quark by 
QGP. Hence the properties of the QGP can be 
extracted from the suppression factor. 



A. Nuclear suppression factor 

The variation of the nuclear suppression fac- 
tor, Raa [3 ■^ith Pt of the electron orig- 
inating from the decays of D and B mesons 
have been displayed in Fig. [7] for RHIC initial 
condition [Ti — 300 MeV). A less suppression 
of B is observed compared to D. The theo- 
retical results show a slight upward trend for 
Pt above 10 GeV both for mesons contain- 
ing charm and bottom quarks. Similar trend 
has recently been experimentally observed for 
light mesons at LHC energy [22] • This may 
originate from the fact that the drag (and 
hence the quenching) for charm and bottom 
quarks are less at higher momentum. 

The same formalism is extended to evalu- 
ate the nuclear suppression factor, Raa both 
for charm and bottom at LHC energy. Re- 
sult has been compared with the recent AL- 
ICE data(Ref. (23,]) in VigM The data is re- 
produced well by assuming formation of QGP 
at an initial temperature ~ 550 MeV after 
Pb-hPb collisions at ,/s^ = 2.76 TeV. 



B. Elliptic flow of heavy flavours 



In Fig. |9] we compare the experimental 



data obtained by the PHENIX [2J] collabo- 
rations for Au + Au minimum bias collisions 



at -y/sNjv = 200 GeV with theoretical results 
obtain in the present work. We observe that 
the data can be reproduced by including both 
radiative and collisional loss with Cg ~ l/\/4. 
In this case V2^ first increases and reaches 
a maximum of about 7% then saturates for 
Pt > 2 GeV. However, with ideal equation of 
state {cl = 1/3) we fail to reproduce data. 
This is because with larger value of c^ the 
system expands faster as a result has shorter 
life time for fixed T,; and Tc- Consequently 
the heavy quarks get lesser time to interact 



1.2 
1 
0.8 
0.6 
0.4 
0.2 




- D 
- B 
• D°-ALICE 




6 8 10 12 14 16 18 20 

Pt (GeV) 



FIG. 8: Raa as a function of pr for D and B 
mesons at LHC. Experimental data taken from 
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FIG. 9: Elliptic flow of single electrons originating 
from the heavy mesons decays. 



with the expanding thermal system and fails 
to generate enough flow. From the energy dis- 
sipation we have evaluated the shear viscosity 
(ry) to entropy (s) density ratio using the rela- 
tion HI: t{/s - 1.25TVoT;vhere q = {pl)/L 
and dE/dx ~ as{px) [221; L is the length 
traversed by the heavy quark. The average 
value oirj/s ~ 0.1-0.2, close to the AdS/CFT 
bound [23. 



IV. SUMMARY 

In this work we have discussed the produc- 
tions of lepton pairs from nuclear collisions 
at relativistic energies and shown that lepton 
pairs can trace the evolution of collectivity of 
the system. The elliptic flow and the nuclear 
suppression factor of the electrons originating 
from the heavy flavour decays have been stud- 
ied by including both the radiative and the 
coUisional processes of energy loss in evalu- 
ating the effective drag and diffusion coeffi- 
cients of the heavy quarks. The results have 
been compared with the available experimen- 
tal data and properties of QGP expected to 
be formed at RHIC collisions have been ex- 
tracted. 
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